Astronomy & Astrophysics manuscript no. 


1st February 2008 


(DOI: will be inserted by hand later) 





Extending the limits of globule detection* 

ISOPHOT Serendipity Survey observations of interstellar clouds 
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Abstract A faint 7i7o = 4 MJysr -1 bipolar globule was discovered with the ISOPHOT 170 /an Serendipity Survey (ISOSS). 
ISOSS J 20246+6541 is a cold (T d » 14.5 K) FIR source without an IRAS pointsource counterpart. In the Digitized Sky Survey 
B band it is seen as a 3'size bipolar nebulosity with an average excess surface brightness of ~ 26mag/D". The CO column 
density distribution determined by multi-isotopic, multi-level CO measurements with the IRAM-30m telescope agrees well 
with the optical appearance. An average hydrogen column density of w 10 21 cm~ 2 was derived from both the FIR and CO 
data. Using a kinematic distance estimate of 400 pc the NLTE modelling of the CO, HCO + , and CS measurements gives a peak 



density of i 



. The multiwavelength data characterise ISOSS 20246+6541 as a representative of a class of globules 



which has not been discovered so far due to their small angular size and low 100/im brightness. A significant overabundance of 
13 CO is found X( 13 CO) > 150 x A"(C 18 0). This is likely due to isotope selective chemical processes. 
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1. Introduction 

Bok globules were originally detected in absorption against 
HII regions by Bok and Reilly (1947), and are known as 
small, dense interstellar clouds in the solar neighbourhod (d iS 
400pc). They were identified optically, and mostly the nearby 
ones have been catalogued so far. It is expected that they are 
similarily common elsewhere in the Galactic disk and in fact 
a few distant ones are already investigated eg. by Launhardt 
& Henning (1997). Their FIR properties were determined by 



Clemens & Barvainis (1988) based on IRAS data. They empha- 



sized the importance of finding FIR faint cold globules since 
these are representatives of the inactive (i.e. non-starforming) 
interstellar medium. Starless globules with small apparent size 
and low temperature (thus also very low 100 /mi brightness) 
can be seen only by good sensitivity at wavelengths over 
100 /im. An ISOPHOT study of pre-stellar cores in dark clouds 
with low 100/im brigh tness were recently reported by W ard- 
Thompson et al. (2002). The ISOPHOT 170 /im Serendipity 
Survey (ISOSS) can be used to locate cold galactic objects even 



without any preliminary identification. This raises the possi- 
bility of detection of the "missing" globules via ISOSS. We 
present our results on ISOSS J 20246+6541 proving that it is, 
indeed, a small and cold isolated molecular cloud - one of the 
so far missed population. 

2. Observation and data analysis 

We searched ISOSS at medium galactic latitudes for 170 /im 
point sources without IRAS counterparts. The ISOSS measure- 
ments, calibr ation, and data ana l ysis for interstellar c louds are 
described by [Bogun et al. (1996| ), |Stickel et al. (2000| ) and Toth 
et al. (2000) respectively. ISOSS J 20246+6541 was detected 
as an ii 70 (peak) = 4MJysr~ 1 (> 4ct) pointlike (FWHM 
< 2.5') source at RA(2000)=20 h 24 m 36 s Dec(2000)=+65°40' 
(1=99.80°, b=15.70°). It is without an IRAS point source coun- 
terpart. 

The optical appearance of CISS1 and its neighourhood 
was studied using 40' x 40' Digital Sky Survey (DSS) blue 
and red images. Photometry on DSS plates of all USNO 
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(Monet et al. 1998) sources have been done where R and B 



< 18mag. Integrated photographic density was derived for all 
these USNO sources. Photographic density values DN of the 
background were measured at each plate at 25 positions where 
neither star, nor significant background emission enhancement 



2 



L.V. Toth et al.: Extending the limits of globule detection 




-4 -2 
v lsr [km/s] 

Figure 1. CO spectra at the 13 CO(1-0) peak, (2-1) lines are 
overlaid as filled histograms, a, 12 CO(1-0) and (2-1) spectra, 
b, 13 CO(1-0) and (2-1) spectra, c, C 18 O(l-0) and (2-1) spec- 
tra. 



was seen. The background value is DN(BG) R = 4590 ± 190 
and DN(BG)b = 4610 ± 230 in photographic density units 
for R and B bands respectively. Photometric calibration of the 
plates was made following Cutri's (1993) method. The photo- 
graphic density to magnitude calibration formulae were derived 
as logio(DN) vs. magnitude relations. The scatter is large for 
faint stars, and the extrapolation to low photographical density 
(i.e. high magnitude) values results in a large error bar. In or- 
der to make surface brightness maps of ISOSS J 20246+6541, 
stars were removed from the optical images by substituting the 
average surrounding DN value. 

The J =(1-0) and (2-1) rotational lines of CO, 13 CO and 
C ls O were measured in Aug. 1998 with the IRAM-30 tele- 
scope. A 4.5' x 4.5' region centered on ISOSS 20246+6541 
was mappen in on-the-fly mode in 12 CO(1-0), 12 CO(2-l) 
and 13 CO(1-0) with angular resolutions of 20", 10", 20" re- 
spectively. Pointed measurements were performed in all the 
above mentioned transitions at the 13 CO(1-0) peak position 
(22 h 24 m 46 s , 65°40'01"). Spectral resolutions of 0.1, 0.05 and 
0.03kms _1 were used. These spectra are shown in Fig. |l[ The 
lines are narrow with a FWHM of < 0.5kms _1 . The result of 
the 1998 C ls O pointed measurements were confirmed in Sept. 
2000 when also the CS(2-1) and (3-2) as well as the HCO+(l- 
0) and (2-1) transitions were observed with the IRAM-30 tele- 
scope. The C ls O, CS(2-1) and HCO+(1-0) lines were de- 
tected with S/N > 5 and well resolved with the 0.03kms _1 
spectral channels. The data were calibrated to the scale of cor- 
rected antenna temperature, T£, by observing loads at ambient 
and cold temperature, as in the conventional "chopper-wheel" 
calibration for millimetre wavelength observations. The data 
were converted to the main-beam scale by applying the relation 



Tmb = ?a x (Bes/Fefi)- The values for the main-beam effi- 
ciency, B cS , are 0.70 and 0.42 (July 1998), and 0.75 and 0.53 
(September 2000), and the values for the forward efficiency, 
F cS , are 0.92 and 0.85 for CO(1-0) and CO(2-l), respectively. 
We checked the calibration by observing the standard sources 
DR 21 and W30H and found it always to be better than 20 %. 

3. Results 

3. 1. FIR results and derived parameters 

ISOSS J 20246+6541 is a faint FIR source with an upper limit 
on the Tioo brightness of about 0.5MJysr _1 derived from the 
"raw" (IRD S format) IRAS data Comparison of ISOSS and 
ISSA data (Wheelock et al. 1994) was made as described by 



Toth et al. (200C ). An upper limit of the colour temperature of 
14.5 K was estimated from the bisector slope of the Iuo vs. 
Jioo scatter plot. Assuming a dust temperature of 14.5 K, an 
average dust column density of 1.3 x 10~ 5 gcm~ 2 was derived 
following Hildebrand (1983). This corresponds to an average 
hydrogen column density of /V(HI+2H 2 ) = 8.6 x 10 20 cm~ 2 
where a hydrogen-to-dust mass ratio of 110 (Launhardt & 
Henning 1997) was assumed. 

3.2. Optical images 

ISOSS J 20246+6541 appears as a faint, isolated reflection 
cloud west of the L1155/L1157 cloud complex. When 
smoothed to 15" the excess diffuse surface brightness distri- 
bution of ISOSS 20246+6541 shows a "bright" lobe at the NE 
and a fainter fragmented one at SE. The POSS B band diffuse 
surface brightness is shown in Fig. ||a, where the contours are 
drawn from 22.7mag/D" by — 0.05mag/D" . The lowest con- 
tour corresponds to three times the standard deviation of pho- 
tographical density outside the globule on the star-removed, 
smoothed image. The average backgr ound surface brightn ess 
is 22.8/D". Linear interpolation of the ^einert et al. (1998 ) ta- 
ble of the B band sky brightness distribution gives at the the 
cloud position 102 Sw, which agrees well with our result. The 
average excess surface brightness of ISOSS J 20246+6541 is 
w 26mag/D"in the DSS B band. The same morphology is 
shown by the B and R surface brightness distributions, except 
that the SE lobe is slightly more red than the NW lobe. The 
lack of a dark core indicates an average visual extinction below 
1.5 mag. 

3.3. Molecular line results, and derived parameters 

All observed lines show a LSR velocity of -2.7kms _1 . The 
FWHM widths of the 12 CO lines are around 0.4kms _1 , and 
around 0.3kms _1 for all the other detected lines. The 12 CO and 
13 CO integrated intensity distributions are shown in Figs. ^ b 
and c. The bipolar shape is well seen and the CO line intensities 
are in accordance with the excess surface brightness distribu- 
tion, the NE lobe being much brighter. The two lobes show the 
same radial velocity. 

The physical parameters were at first approximated from 
the 12 CO(l - 0) and 13 CO(l - 0) spectra of the peak in- 
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Figure 2. Dust and gas in ISOSS 20246+6541, reference coordinate (0,0): RA=20 h 24 m 44 s DEC=65°40'04" 

a: B band diffuse surface brightness, contours are from 22.7mag/D" (i.e. background level + 3 x a converted to magnitudes) 

-0.05mag/D" 

The diffuse surface brightness contours indicate a sharp cloud edge at NE. The structures at offset coordinates (-0.5, -1.8) are 
less reliable due to residuals of a removed 11.9mag star GSC 04245-00803. b: W( 12 CO (1-0)), contours are from 1.5Kkms _1 
with 1.5 Kkms -1 steps; c: 13 CO (1-0), contours are from l.OKkms -1 with l.OKkms -1 steps; CO isotopomer line intensities are 
integrated in the velocity interval -5 kms -1 < vlsr < Okms -1 . 



tensity positions, assuming an isothermal spherical cloud in 
LTE which uniformly fills the beam. The LTE CO(1-0) exci- 
tation temperature is T cx (CO) = 9.2 K, characteristic of non- 
starforming molecular clouds. We derived an optical depth of 
r( 13 CO) = 0.87 from the observed CO(1-0) and 13 CO(1-0) 
line ratio assuming the terrestrial isotopic density ratio of 89. 
The CO column density was calculated using the standard LTE 
formula: 

iV( 13 CO) = 2.57 x 10 14 T kin i T(13( / 0) t v „ resulting in 

CXPl > fcT kin I 

iV( 13 CO) = 8.9 x 10 14 cm~ 2 , which corresponds to iV(H 2 ) rs 
10 21 cm- 2 assuming N(U 2 ) « 10 6 iV 13 (CO). This result is in 
agreement with the column density derived from the FIR data. 
The 12 CO and 13 CO lines trace the ISM well over most of the 
cloud. A comparison of C ls O to 13 CO lines (both 1-0 and 2- 
1) at the centre of the NE lobe indicates an underabundance of 
C ls O by a factor of 4. This effect is expected in cold clouds 
with moderate density, exposed to UV radiation (Glassgold 
et al. 1985). A more careful modelling may account for it as 
we show in the Discussion session. 



4. Discussion 

4. 1. Distance of ISOSS J 20246+6541 

The distance of ISOSS 20246+6541 can be estimated, relat- 
ing it to its neighbours. Its nearest neighbours are LI 122 and 
(|Lynds 1998|), and the YDM97 COl (|Yonekura 1997|). The 



CO survey of Cepheus by Yonekura (1997) covers the po 



sition of ISOSS J 20246+6541 and their Fig. 6 a indicates few 
small clouds around ISOSS J 20246+6541 (i.e. YDM97 COl, 
YDM97 C02, YDM97 C03) and one even smaller unnum- 
bered peak very close to ISOSS J 20246+6541 at 1= 100? 1 b= 
15? 7. All the listed Yonekura clouds have Ulsr > +2.5kms~ 1 , 
and they are counted into the "close group" of clouds, which on 



the other hand is associated with extended FIR features around 
ISOSS J 20246+6541. The nearest molecular clouds with neg- 
ative ulsr are YDM97 C07, C09, CO10 at I w 103?0 b w 
16?7. ISOSS 20246+6541 itself has u LSR = -2.7kms- 1 . It 
probably belongs to the ISM layer of the nearby Cepheus Flare 
GMC, and is located at about 400pc ( |Kun 1998|). We note that 
applying the size-linewidth relation of Larson (1981) the glob- 
ule may be between 100 and 400pc. 

4.2. Radiative transfer models 

We have modelled the NE lobe of the bipolar globule with 
spherically symmetric cloud models, although the NE clump 
shows some deviations from spherical symmetry in both 
12 CO and 13 CO (see Fig. | b and c. With RA=20 h 24 m 44 s 
Dec=+65°40'04" as the centre position, we have averaged 
spectra in concentric rings with radii increasing by 10" inter- 
vals up to a radii of 90". The effective resolution of the av- 
eraged spectra is 40" for the J =1-0 lines and 20" for the 
J =2-1 lines. 

We set the cloud parameters as follows. 

(1) We assume a density distribution n ~ r 
ratio 100 between the centre and the cloud surface. 

(2) The kinetic temperature is assumed to rise linearly from 
the cloud centre. This is a crude approximation of the actual 
temperature structure of for a small, spherically symmetric 
globule without internal heating sources (e.g. Leung 1985; 
Nelson & Langer 1999) but will suffice for the present pur- 
poses. The temperature gradient, i.e. the difference between 
the outermost and innermost shells AT = OK, 6 K or 10K. 
Higher contrast means too high a temperature for the outer 
cloud, in contradiction with the observed small linewidth. 

(3) Extinction-dependent relative molecular abundances 
X (molecule) = "( I ° ol g :ulc ) were estimated according to 



-1.5 



with a density 
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Warin et al. (1996 ). The cloud is cold, exposed to UV radiation 
and it has a peak visual extinction between 1 and 2 mag. In 
these conditions isotope selective processes result in a relative 
overabundance of 13 CO and relative underabundance of C ls O 
according to Bally and Langer (1982). When applying the 
Warin et al. (1996) relative abundances, we introduced an 
intrinsic extinction at the cloud boundary since the 12 CO 
lines are not vanishing at the boundary of the NE lobe. This 
assumption is supported by the presence of surrounding 
extended cirrus-like emission seen in 100 /im on the ISSA 
image. 

(4). Distance: 100, 200, 400, 600, 800, 1000, and 2000pc were 
tested. 

When the density, temperature, relative abundance distri- 
butions and the distance is set to a value allowed by the above 
constrains, the free model parameters are the central density 
(n c ), the intrinsic linewidth (Av) and the angular diameter (D) 
of the model cloud. The radiat ive transfer pr oblem is solved 
with Monte Carlo simulation (luvela 1997). The computed 



spectra are convolved to the resolution of the observed spec- 
tra and the quality of the fit between the two is estimated with 
a weighted \ 2 value. The model cloud is divided into 3 1 shells 
of equal thickness and the free parameters are optimized sep- 
arately for 12 CO and 13 CO. We then select the set of parame- 
ters which provides the best fits for both. Since the 12 CO ob- 
servations only probe the outer layers of the cloud, the 12 CO 
based estimate of the column density is uncertain. Modelling 
based on the 13 CO line, however, gives surprisingly similar 
results when the appropriate average relative abundance value 
X( 13 CO) ~ 1.1 x 10 -6 is selected. The models are not sen- 
sitive to 20% changes in the average molecular abundances or 
density, although similar changes of the kinetic temperature or 
size are, critical (see Fig. ph. 
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from the 12 CO 13 CO analysis. The relative abundances were 
varied up to 100% and the other parameters up to 30%. The 
C ls O lines were best reproduced assuming an average rela- 
tive abundance of X( 13 CO)/X(C 18 0) « 150, an extreme 
but possible underabundance by factor of 28 (Glassgold et al. 
1985). The pointed measurements supported the density and 
temperature results shown in Fig. ||. The derived NLTE ki- 
netic temperature is around UK and assuming a distance of 
400 pc the peak hydrogen density and the size of the NE lobe is 
n c =6.7-10 4 cm -3 and 0.12pc respectively. The column density 



estimate is A^H^) 
be w 3M©. 



2 • 10 cm . The total gas mass would 



With 7V(H2) ~10 cm , the visual extinction towards 



the cloud centre is A v «1 (Bohlin et al. 1978) and for UV 



Figure 3. Results of the radiative transfer modelling of the 
northern lobe assuming distances: lOOpc (yellow), 200pc 
(lilac), 400pc (red), lOOOpc (blue). The relative \ 2 minimum 
regions are plotted as shaded regions in selected planes of the 
parameter space. Two further contours of the relative \ 2 of 1 . 14 
& 1 .29 are overlaid. (Extinction dependent relative molecular 
abundances were used.) a: density versus kinetic temperature 
(left); b:density versus size (right) 

Synthetic spectra for C ls O CS and HCO + were generated 
with the NLTE model using the best parameter sets (lowest x 2 ) 



photons the cloud is optically thick unless it is very clumpy. 
External heating, however, is reduced by the surrounding ISM, 
which is represented by the nonvanishing 12 CO lines. This may 
be the reason that a moderate 6K temperature contrast was 
found to be more likely than a 10K contrast or isothermal. 

Although the SW clump is clearly elongated similar mod- 
elling was carried out for that part of the cloud. Observed spec- 
tra were averaged over annuli at radii up to 50" from the clump 
centre. Assuming a model where the kinetic temperature in- 
creases linearly from the centre, we obtain a peak column den- 
sity of A f (H 2 )=2.9-10 20 cm" 2 based on the CO spectra. 

5. Concluding remarks 

(i) ISOSS 20246+6541 is a small isolated molecular cloud dis- 
covered by its 170 fim emission. 

(ii) The globule appears as bipolar in both optical surface 
brightness and mm-line maps. 

(iii) A distance of 400 pc is likely, setting the diameter of the 
northern lobe to approximately 0. 1 pc. 

(iv) The average physical parameters are: Tki n = UK, 
n(H 2 ) = 7000cm- 3 . 

(v) In the 5° environment of ISOSS 20246+6541 we found fur- 
ther 15 ISOSS sources with similar FIR parameters. One 
of those is another previously unknown globule as seen on 
DSS2 images. 

(vi) The population of small and faint starless globules can only 
be explored by high sensitivity FIR measurements such as 
those carried out by PHT-C2 on board ISO. 

(vii) We note that there are no other ISO measurements but the 
ISOSS within 40'distance from ISOSS J 20246+6541. Our 
finding thus also demonstrates the importance of such ubi- 
ased surveys by future missions like Planck. 
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